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Abstract 

After inflation occurs, a weakly coupled scalar field will in general not be in 
thermal equilibrium but have a distribution of values determined by the in- 
flationary Hubble parameter. If such a field subsequently undergoes discrete 
symmetry breaking, then the different degenerate vacua may not be equally 
populated so the domain walls which form will be 'biased' and the wall net- 
work will subsequently collapse. Thus the cosmological domain wall problem 
may be solved for sufficiently weakly coupled fields in a post-inflationary uni- 
verse. We quantify the criteria for determining whether this does happen, 
using a Higgs-like potential with a spontaneously broken Z2 symmetry. 

11.27.-Fd, 98.80.Cq, 05.10.Gg 



1 



Typeset using REVT^ 



I. INTRODUCTION 



It is generally believed that if discrete symmetries of scalar fields are spontaneously broken as 
the universe cools down, then there would be severe difficulties for its subsequent evolution [||. 
This is because topological defects — domain walls — would form at the boundaries of the different 
degenerate vacua chosen in causally disconnected regions following the symmetry breaking phase 
transition and would eventually come to dominate the total energy density, in conflict with 
observations . To avoid this requires the energy scale of the symmetry breaking phase transition 
to be lower than ~ 100 MeV; in fact it must be less than ~ 1 MeV if the anisotropy induced by 
the walls in the cosmic microwave background radiation is to be below experimental limits 
This is a severe constraint on attempts to extend physics beyond the Standard Model, which often 
involve introducing such discrete symmetries [Q]. 

In some special circumstances, the broken discrete symmetry may not be restored at high 
temperature so domain walls would never form ||5|. There would also appear to be no problem if 
the symmetry breaking occurs prior to inflation since one would then expect the density of any 
resulting topological defects to be exponentially diluted away. However defects can still form in 
this case through quantum fluctuations of the scalar field during inflation ||^-^ if the mass of 
the field is less than the inflationary Hubble parameter. When inflation is driven by a F-term 
supergravity potential, the breaking of supersymmetry by the large vacuum energy gives all scalar 
fields, including the inflaton itself, a mass-squared of 0{H^) ||l^,14|; in this case fluctuations are 
negligible and walls will not form. However if we consider instead e.g. a D-term or no-scale 
inflationary potential j^, the scalar field may remain light relative to the Hubble parameter so the 
above mechanism will be operative and domain walls will form. 

Although the first such paper |^] considered an axion field, subsequent work [^J9| has been 
mainly concerned with scalar fields which have sufficiently strong couplings that the vacuum ex- 
pectation value (vev) during inflation does not increase much above the Hubble parameter. The 
field then remains uncorrelated on spatial scales larger than the Hubble radius and defects form 
during or at the end of inflation. However very weakly coupled scalar fields are arguably of more 
interest in cosmology. For example the field responsible for driving inflation should have very small 
couplings in order that its quantum fluctuations not contribute excessively to the anisotropy of the 
cosmic microwave background [p. There has been much interest in 'quintessence' |10| — a very 
weakly coupled evolving field that may account for the tiny vacuum energy that is suggested by 
some astronomical observations. Weakly coupled fields can also be a source of dark matter through 
their coherent oscillations [11|. In a recent paper |12| an extremely weakly coupled dilaton field 
that forms domain walls is proposed as a way of binding the matter in spiral galaxies and producing 
their characteristic flat rotation curves (as an alternative to cold dark matter). Particularly in the 
context of this model, it is interesting to ask whether the above mechanism would indeed create 
stable domain walls. 

The point is that such a very weakly coupled field will be correlated on super-horizon scales 
at the end of inflation and not be brought back into thermal equilibrium during the reheating 
process since it has no couplings to the thermal plasma or to the inflaton. The field will oscillate 
coherently during the post-inflationary Friedman-Lemaitre-Robertson- Walker (FLRW) expansion 
era and when the expansion redshift reduces the energy in its coherent oscillations it will settle into 
different symmetry-breaking vacua on spatial scales larger than the Hubble radius, thus forming 
defects. It would be likely for the same vacuum to be chosen in different (apparently causally 
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disconnected) regions. A 'bias' could thus be generated in the probabihties for populating the 
distinct vacuua even if they are energetically degenerate After the walls form, such a bias, 

even if very small, will result in exponential decay of the wall network, as has been demonstrated 
both analytically and numerically |16-|l^. Thus there may be no domain wall problem for weakly 



coupled fields in a post-inflationary universe. 

Our aim is to quantify the bias that would be created for such a hypothetical field with specified 
properties in order to determine the fate of the domain walls formed. We consider the problem 
in its simplest form and focus on domain wall formation through inflationary fluctuations in a 
spontaneously broken theory of a real scalar field with the Higgs-like potential 

V[4>) = -\n^^^\M' = \[f-v^)\ (1) 

where v = m/y/X. First we study (Section ||) the stochastic evolution of the field perturbations 
during the inflationary epoch which set the relevant initial conditions. In Section Jl^ we follow the 
evolution of the fluctuations during the FLRW expansion until the field drops into its potential 
minima and the domain walls are formed. The bias between the degenerate vacua is calculated in 
Section [V . We review the history of the field evolution in Section ^ and identify the regions in 



the parameter space of the above model where domain walls do not survive. Finally we present 
our conclusions and comment on specific models concerning domain walls such as Ref. [12|. 



II. STOCHASTIC APPROACH FOR THE INITIAL CONDITIONS 

During inflation the smooth component of a slowly evolving scalar field can be considered (on 
scales larger than the horizon) to be a classical variable subject to stochastic noise (contributed 
by the field modes whose exponentially increasing wavelength causes them to 'exit the horizon', 
becoming part of the coarse-grained field) p9|- p2| , p^ . The Langevin equation governing the coarse- 
grained field (j) is 



where the white noise r] satisfies 

{r^{t)r^{t')) = 6{t - t') . (3) 

This equation can be restated as a Fokker-Plank equation for the probability distribution P{(j), t) 
of the field values in a given coarse-grained domain |jT9|: 

-P(0,t) . (4) 



dt 5(/> \3Hi d(j) V 87r2 

Here H[, the Hubble parameter during inflation, is taken to be independent of (/>, i.e this field is 
assumed not to contribute significantly to the vacuum energy during inflation. (In the analogous 
equation for the inflaton field, Hi is itself a function of (p giving rise to ordering ambiguities in the 
corresponding Fokker-Planck equation.) 

If the field evolution is 'slow-roll', then the force term (involving the potential) can be neglected 
compared with the noise term in Eq. (^) . Then starting from a given value of the field (j) = (p averaged 
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over a patch of size ^ when cosmologicahy relevant scales 'exit the horizon', the solution to this 
equation is the Gaussian distribution: 



P^ = P( 



'87r3 

Wt 



exp 



Hh 



(5) 



where Hi is taken to be approximately constant as is required for successful inflation Q. 



The mean value (0^) = H^t/Air^ grows linearly with time as in Brownian motion |25]. Relating 
the time during inflation to the cosmological scale through l^^ ~ H\e~^^*, we can write the 
probability distribution as [^,15| 



27r(T 



exp 



1 

2^ 



(6) 



where o"^ is the quadratic dispersion of the field. This can be obtained, with reference to the 
noise term in Eq.(^), as the sum of independent Gaussian distributions with dispersion-squared 
(i?i/27r)^, one for each e-fold of inflation. The sum is to be taken over the period when scales 
between /min and /max leave the horizon, where l^m corresponds to the Hubble radius at some 
moment during the FLRW era (representing an ultraviolet cutoff, given that we are interested in 
the super-horizon behavior) and /max 

is the biggest spatial scale of interest, i.e. of order the present 
Hubble radius H^"^ . The total dispersion-squared is then just the sum of the dispersion-squared 
for the independent probabilities: 
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(7) 



i.e. of 0{Hf ) in the cases of interest. 

The formula assumes that the value of the force term in Eq.(^) is negligible in comparison 
with the noise term so the value of (j) is not determined. However, if inflation continues for a large 
number of e-folds the force term will impede the tendency of the distribution to widen indefinitely. 
In this case stochastic equilibrium is achieved and it is possible to give a probabilistic prediction 
for the initial (f) using the stationary solution for the Fokker-Planck equation (while Eq.(|6|) still 
gives the distribution for the field at the end of inflation on cosmologicahy interesting scales). Q 

The stationary case dP/dt = can be solved to obtain the probability distribution of the 
averaged field (p: 



Ci exp 



Svr- 



3 i 



-I- C2 exp 



8tt' 



exp 




(8) 



If the potential is an even function the first term is also even while the second term is odd. 
Therefore, if the potential remains positive for large (j), the second term will have greater absolute 
values than the first at some point, and as it is an odd function the probability would have negative 
values. This shows that the second term is unphysical. Thus, in the stationary case the normalized 
probability distribution is just: 



^This has been shown in studies of 'eternal inflation' |26|. 
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Pi = exp 



Stt^ V{(t)) 
3 Ht 




(9) 



Defining the dimensionless normalized fields x = 0/^^) ^"^^ X = 4'/'^, this is a one-parameter 
function for our chosen potential (||): 



4' 3^ 



exp 



, f3 = \ 



^1/4 -^i 



(10) 



where /(x, y) is the modified Bessel function of first kind. This is plotted in Fig.|^ for various values 
of /3. 
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FIG. 1. Stochastic equilibrium probability distribution for the (normalized) field vev during 
inflation. The curves shown are for (3 = 1.2, 1.6, 2.8, 5 in order of decreasing height. 

Since the height of the potential barrier at the origin separating the two vacua is = \\v^ 
from Eq.(||), the parameter (3 in Eq.(lO) is just /3 = i/i/\/2/i. Thus for small (3 the probability 
of transitions between the potential wells is exponentially suppressed; the discrete symmetry is 
broken during inflation and the field is localized in one minimum or the other so domain walls will 
not form (assuming that the symmetry is not restored during reheating after inflation). However 
when H[ ^ h the fluctuations are large enough to make the probability distribution flat as seen 
in Fig.|l[ 

In general the average value of the potential energy in the field (p is 



(11) 



{V)= Vi<P)P^d4 = giP)Ht 



where g{P) is a smooth positive function with maximum value 0.02 and limiting values 1.9 x 10 ^ for 
/3 and 9.4 x lO'^ for P ^ oo. Thus the energy density of the scalar field under consideration is 
indeed negligible when compared with that of the inflaton Vi ~ 3H?Mp, as was implicitly assumed 
in neglecting the (f) dependence of Hi. (Here Mp = l/\/87r G ~ 2.4 x 10^^ GeV is the normalized 
Planck scale.) Note that the COBE observations of large-scale anisotropy in the cosmic microwave 
background set a strict upper bound on the Hubble parameter during inflation |0] 
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V^'^ « 0.027Mp ^ ^ < 4.2 X 10^^ (12) 

This still allows the energy scale of inflation to be as high as the GUT scale (~ 10^^ GeV) but in 
specific models it can be much lower than this, in particular in 'new' inflationary models with a 
quadratic leading term in the potential |28|. 



The average value of is just v for small /?, while for large /? it is ~ 0.305 Fi/A^/^. Thus if 
A ^ 8 X 10^^, the field vev grows above H\. The number of e-folds of inflation that are necessary 
to achieve the stationary distribution in this case must exceed {{(p) /6(f))'^ ~ A~^/^, taking the field 
increment per e-fold to be Sep — Hi/In. Otherwise one cannot predict an unique probability 
distribution for the field. 

The slow-roll condition (j) <^ 3Hi(/) implies <C 9H^. For the potential (|l|) and the 

distribution ( [Tol) this translates into <C f-fff if the quadratic term in the potential dominates 
(small (3). When the quartic term dominates (large (3) one must require A ^ 1. If the slow-roll 
conditions are not obeyed, the classical stochastic treatment does not apply and the quantum 
creation of particles is exponentially suppressed |2^. For the case under consideration, h <^ Hi, 
the relevant slow-roll condition is A ^ 1 and these two conditions together imply m ^ H\. 



III. THE EVOLUTION DURING THE FLRW ERA 

In Ref. [15 1 it was shown that the probability distribution for the coarse-grained field outside 



the horizon during the FLRW era following the infiationary era evolves into a non-Gaussian distri- 
bution. We are interested here in obtaining a discrete probability distribution for the two vacua of 
the field. We follow the evolution of the field from its initial value (po at the end of infiation until 
it settles down into one of its discrete minima. By obtaining the distribution function f{(j)o) whose 
values are ±1 according to the vacuum finally chosen, we can compute the bias between the vacua. 

The equation of motion for the modes of a scalar field outside the horizon during the FLRW 
era is 

d'^d) c dd> , 

+ 3--- + y' = 0, where c = Ht . (13) 
dt^ t at 

We have rewritten the Hubble parameter in terms of the variable c which equals 1/2(2/3) for a 
radiation- (matter-) dominated universe. The initial conditions are c/to = Hi, (p = (po, and (pQ = 0. 
In general there is no exact solution for this equation, however a good approximation may easily 
be obtained [p7| ]. 

Let us first suppose that the friction term is relatively unimportant so at first approximation 
the field oscillates in its potential minimum at the end of infiation according to the Lagrangian 
L = ^(j)'^ — V. Then energy conservation implies 



/,2 



2(l^max - V) , (14) 



where Vmax is the maximum value of the potential energy during the oscillation. Thus the oscillation 
period is given (for a symmetric potential) by 

At = 4 / , ^ . (15) 

Jo \/2(Knax - V) 
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When H <^ Lo = 2-7r/At the friction term in Eq.(p!^) is lo/H times smaller than the other terms so 
the assumption of negligible friction is valid. Until H drops below uj the field remains approximately 
fixed since the friction is relatively high and the dynamical time exceeds the expansion age H~^. 
After this point the field starts oscillating, losing energy density of order 

3H<j)d(t> = UH p{V^^^ - y)d0 (16) 

in each oscillation. Here we have used the fact that in this regime H is approximately constant 
during each oscillation period. 

For future use we also consider the case of a general power law potential 

V{<t>) = -<t>^. (17) 
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Let (/>n be the values of 4> that are turning points of the trajectory, tn the corresponding times and 
Pn the corresponding energy densities. Evidently and tn grows, eventually going to infinity. 

Equations (|l^) and (|l6|) imply the following relations between these quantities 

= = (18) 

Atn = fci 01-^/2 , (19) 

where h = (^^r^^+TT^f^) -^"^^^ ^2 = c ^(2^-)!^('i/2+i/7) ) A"^/^. The solution to these 
recurrence relations are the power-laws: 

6c 

'^" = ^^""' " = - 2(l + 3c)+7(l-3c) ' 
.b u (2 + 7) 



tin, b 



2(1 + 3c) +7(1 -3c) 



The energy density is ^ ^ t'^'^/^ . In terms of the cosmological scale- factor R this can be 
written, p^ ~ R^'^^''^, i.e. 

p^ ~ i?-67/(2+7) ^ (21) 



independently of c |27]. The number of oscillations goes as n ~ (j)^^"" ~ 

t^/b ^ iji/cfe. Table pi 

shows the exponents of different quantities as functions of n, R and t for the cases of radiation- 
and matter- domination and for quadratic and quartic potentials. 

Let us return to the double well potential (|^) which interests us here. As we saw in the preceding 
Section, the quartic term dominates at the end of inflation in the cases of interest, i.e /3 3> 1. The 
subsequent evolution starting with (p = (pQ at t = to goes as follows. The field oscillates in the 
quartic term-dominated potential (where the mass term can be neglected) , decreasing in amplitude 
due to friction caused by the Hubble expansion, until V^^^ drops below the height of the barrier h 
separating the two vacua, or in other words 4> becomes of 0{v), and the field settles down in one of 
its potential minima. The energy density in oscillations then decreases as for radiation {p^ ~ R~^) 
independently of the rate of expansion. The relation between the number of half-oscillations and 
the field amplitude is n ^ (p^^ if the universe is radiation-dominated and n ~ (p^^^"^ if it is matter- 
dominated. Thus (p will be found around one of its vevs ± v after having completed n = t^Xo 



7 



half-oscillations in the radiation-dominated case and n = Tm-^/xo in the matter-dominated case, 
where xo = ^o/v, and Tr, Tm are constants. (The numerical values of these constants cannot be 
calculated given the approximations we have made since they depend on the details of the evolution 
towards the end of the oscillations (when the mass term begins to be significant), and even more 
importantly on to (= c/H[), which sets the Hubble parameter at the beginning of the oscillations 
and thus determines the friction term.) In general n is a function of xo; ^0 and A. However, 
rewriting the equation (O) using (f) = and i = \/^(pot, we have: 



dt2 t dt V Xo 



+ h^+{^'-^^)=^- (22) 



This equation has one parameter and the initial conditions are (po = 1, d(/>/dt = 0. Therefore the 
number of oscillations must depend on just xo and on to = c^/X(f)o/Hi. However, as we stated 
before, if initially uj = O.84-v/A0o the field remains frozen. It only begins to oscillate when H 

decreases below w, i.e. at to ~ 1, hence the problem does not depend on to or any other parameter 
as long as we are interested only in the final state. We have checked numerically that n does not 
depend significantly on to when to ^ 0.3. In the present case we have that the initial distribution 
of (1)0 is most probably concentrated around (/jq 0.1 X^^^^Hi, and to ^ 0.1c A~^/'^, so there can be 
no significant dependence on to. (It may be that the initial value of the Hubble parameter for the 
FLRW era. Ho, is somewhat less than its inflationary value, H[, but even in this case the results 
apply for small A.) The proportionality constants can be calculated numerically in this regime and 
are found to be = 0.31 and = 0.63. When to exceeds unity the Hubble parameter is smaller 
than uj at the beginning of the oscillations and the reduction of the friction implies a proportional 
increase of the number of oscillations: Tj-, oc to for to > 1. 

The function f{xo) that equals ±1 according to the value (p = zizv of the final state will then 
change sign just once as n increases by unity. Therefore we can write 

/(xo) = (-1)^"*("'^"\ for c=i, 

^ (_iynt(r^V»^), for c=-, (23) 

3 

where int(a;) is the integer part of x. Small deviations from these expressions occur for low values 
of Xo- The oscillations start when t = t^ 1/\/A</'o — \~^^^H^^ , and whether the universe is 
radiation- or matter- dominated at this time is determined by whether t* is smaller or larger than 
the epoch of matter and radiation equality tgq. (There is also the possibility of a different expansion 
rate during the reheating process.) 
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TABLE I. Power-law exponents for the evolution of field variables during oscillations after 
inflation, for a quadratic (7 = 2) and quartic (7 = 4) potential, assuming a radiation-dominated 
(c = 1/2) and matter-dominated (c = 2/3) universe. 
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Note that the functions (p^) apply only when the oscillations begin and end during a period 
of expansion while c is constant. As an example of a more complex situation consider the case 
where the oscillations start in the radiation-dominated era and end in the matter-dominated one. 
According to Table | the amplitude of the oscillations goes as (/> 1/R. The amplitude at the 
time of matter-radiation equality is then 0eq = 4>o{t*/tcq)^^'^ ~ {4>o/tcq\'^)^^'^ , while the condition 
that the oscillations end in the matter-dominated period is (peq > v. The number of oscillations 
in the radiation-dominated period is just Uj- ~ 4'o/(pcq = (^oqV^'/'o)^''^- Therefore the number of 
half-oscillations completed in the radiation-dominated period is proportional to -y/^o ■ With respect 
to the remaining oscillations that occur in the matter-dominated period, since they start at t^q 
well inside the low friction regime, the constant in Eq.(p3|) has to be scaled by tgq = V^'Pcqteq- 
Accordingly, the matter-dominated period gives the number of oscillations nm — V^fpeqieq^/Xeq = 
Thus, for such mixed situations, different functional dependencies of (/)o are 
expected in the formula for the number of oscillations. 



IV. BIAS 

The bias, defined as the difference in the probabilities of populating the two discrete vacua, is 
given by the convolution 



b{x)= J fiXo)PxiXo,x)dxo. (24) 

The probability distribution of the field values at the end of inflation P^{xoi X) is given by Eq.@ 
(rewritten in terms of the normalized fields x = 0/f and x = 4>/v), while the function /(xo) that 
gives the sign of the field in the final vacuum state was obtained in the previous Section for different 
cosmological situations. (Note that Pxixo^x) is almost independent of spatial scale in the FLRW 
era since observable scales correspond to only a few e- foldings during inflation.) In the following 
we do a detailed analysis of the bias function for the cases of radiation- and matter-dominated 
universes and then present an approximation suitable for more complicated situations. 

The bias in the observable universe is well deflned by Eq.(p^ but since only the probabilistic 
distribution (|T^ is available for x, the predictions are also in terms of a probability distribution 
for the bias which satisfies 

n = E^x§, (25) 

where x is understood as a function of b (inverse of the function (p^)), and the sum is over the 
different branches in the solution of the equation b = 6(x)- The cumulative probability for the bias 
to be less than some particular value is then just the integral: 

P{\b\<x) = Y, P^dx. (26) 

The problem has basically two parameters, 

a = ^~^, (27) 

V V 

which gives the width of the Gaussian probability distribution (^) in terms of the normalized field 
Xo = 0o/^) ^-nd P ~ A~^/^a (using Eq.(lO)) which gives the width of the distribution (10) of the 
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normalized average field x = (p/v- The conditions we are assuming for slow-roll m ^ H[ and A < 1 
imply that a < min(/3, but they do not constrain this parameter to be greater or smaller than 
unity since we also require /3 ^ 1 (i.e. H[ ^ h) in order for the field to be able to jump the 
potential barrier during inflation. 



1. Radiation dominated universe 

In this case the number of half-oscillations n ~ TrXO; and f{xo) = (— ^ square- 
wave. It is easy to see that the bias is then a periodic function of x with period 2t~^, and can be 
expanded in the Fourier series 



4 sin[(2n + l)7rrryl 
- , -,^ exp 



vr 



n=0 



(2n + 1) 



-(2n + l)^a 



(28) 



which converges exponentially fast. For a ^ 1 (i.e. Hi ^ v) we can approximate the series by its 
first term 



b{x) = ^ sin(7rTrx) exp (^--^^ 



(29) 



showing that the bias is a sine function with an exponentially damped amplitude. In the opposite 
case a ^ 1 (i.e. Hi <^ v), more terms of the series must be added so it converges to the Fourier 
series for the square wave, i.e. identical to Eq.(|2^) without the exponential factor. In this limit 
the Gaussian distribution Pxixo^x) is appropriate inside the regions where /(xo) has one sign or 
the other, but moving x from one of these regions to another where f{xo) has opposite sign causes 
the function b{x) to step as the error function, erf(x) = Jq dte~^ . Therefore for q <C 1 we 



can write the periodic function b in the interval [ 



Hx) 



erf ( — ^ 



erf 



1 ix + rr') 

V2 a 



erf 



1 (x 



V2 a 



(30) 



This function is very close to +1 or -1 except in the neighborhood of the origin (or the points nr^ ^) 
where it is linearly dependent on x- 



Kx) 



2X 



vr a 



X<a. 



(31) 



In Fig. 0(a) we show the bias function for several values of a. 

When /? ^ 1, we have that (/iq ~ A~^/^i7i ^ v, and the initial values of x will be distributed 
with equal probability in the interval [— TjT^, r^^^]. Thus, in this case there is no dependence on 
f3, and the results are relatively insensitive to the initial probability distribution for X; even if 
stochastic equilibrium is not achieved during inflation. Therefore, the bias probability function 
P{b) can be simply calculated using Pb{b) = with x in the interval [0, r~^/2], yielding 



vr 



4\ ^ 

— ) exp (—TTT^a) — b 
■kJ 

\J^Tra for a <^l, |6| < 1 



4 / TTT 

for a > 1 , 161 < — exp -^a 

vr V V2 



(32) 
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(a) (b) 

FIG. 2. The bias function for different values of a in (a) the radiation-dominated and (b) the 
matter-dominated case. In the left panel, the curves correspond, from top to bottom, to the values 
a = 0.05, 0.1, 0.2, 0.4, 0.6, 0.8, while on the right, the bias is shown for q = 100. 



where we have shown the maximum value the bias can reach (see Fig.^(a)). The cumulative 
probability for the bias is P{\h\ < x) = 2rrx(x), where the function xi"^^ is just the inverse of the 
bias function in the interval [0, r~^/2]. This gives the approximate answer 



P(|6| < x) = — arcsin 
vr 



pvn 



exp 



a 



\V2 

'2TiT^a X for x ^ 1, a <^ 1 . 
This is plotted in Fig. |^(a) for various values of a. 



vr 



for X < — exp 



V2 



a 



>1, 



(33) 



2. Matter dominated universe 



For the matter-dominated universe we have n = t^^/xo, and the bias will not be a periodic 
function of x- Consequently the bias probability function is more sensitive to the initial probability 
distribution of X; in contrast to the radiation-dominated case. Taking the stochastic distribution 
as the initial probability distribution for x the problem has two parameters rather than just one 
as in the radiation-dominated case. 

Of course, the bias function b{x) depends only on a. If a ^ 1 then the bias will be concentrated 
around ±1 since the width of the distribution will not allow different final vacuum states to be 
reached. The function will be ±1 except near the transition points Xn = Tm^n'^, where it jumps as 
±erf(^^^). In fact this will also be the case for values of x greater than (with a ^ 1), because 
at this point the spacing of xo which results in different final states is > a. In the opposite case, 
X ^ a^, significant cancellation take place and the bias will be reduced. An excellent approximation 
in this regime (if a ^ 1) is given by the formula (^), if we allow for variations in the period and 
the amplitude of the sine function according to the effective variation of a with x (taking into 
account the increasing distance between the Xn with n in the matter dominated case), i.e. 
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(a) (b) 

FIG. 3. The cumulative probability in (a) the radiation-dominated and (b) the mat- 
ter-dominated case. In the left panel, the curves correspond, from bottom to top, to the values 
a. = 0.05,0.1,0.2,0.4,0.6,0.8 (as in Fig. ^(a)), while in the right panel the curves correspond, from 
bottom to top, to the values a/y/P = 0.2, 0.45, 0.7, 1.0, 1.6, 2.2, 3.2. 



Kx) = -sgn(x) exp 



vr 



1 

\x\ 



2^2 / 



sm vrr, 



for X ~ 



a 



(34) 



In Fig. §(b) we show an example of the bias function for the matter-dominated case. 

For the case of interest, i.e. (3^1, the mass term can be neglected in Eq.(lO) and the initial 
probability for x writes 



23/431/4^(5/4)^ 



exp 



(35) 



This function is practically constant for x < /5/2 and then falls until x = where it is nearly zero. 
Thus for /? ^ the bias is ±1 with high probability because x ~ /9- In this case we recover the 
linear behavior (33) of the cumulative probability for the bias in the radiation-dominated case, but 
now the spacing depends on (5: 



P{\b\ <x) = V27r 



a 



X , for a: < 1, a <C , (36) 

0.71, as can be calculated 



where the precise coefficient in the identification Tj. ^ ~ y/j3 is '^Je^^l^^j^^ 
from the average value of the inverse spacing. 

For a < P ^ the bias has higher probability between —1 and 1. In this case the probability 



P(|6| < x) is dominated by the sector in the bias function ( p^ ) where x > - exp 



and thus \x\ <Xx = -(^) V log(f x). Then 



1 f TrrmQ V 
m{2V2 ) 



P{\b\<x) = 2 P^dx = l 



for ^/^< 



a , 



(37) 
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where the function T{x,y) = t^~^e~*dt is the incomplete gamma function. The behavior of 
the cumulative probability function changes from linear to the form (^) as f3/a^ becomes smaller 
and in both limits it depends only on this combination of the parameters as shown in Fig. ^(b). 

Let us recapitulate the main features. In the radiation-dominated case there is no dependence 
on P, and the probability of a small bias b is just ~ ba for a ~ Hi/v <^ 1, while for a ^ 1 it 

2 

increases exponentially towards unity as ~ be" . For the matter-dominated universe, the parameter 
that regulates the bias probability and plays the same role as a in the radiation-dominated case 
is a/\f0 ~ X^/^Hi/v. If /3 = \~^/'^H^/v ^ the probability of a small bias is also small, 

~ ba/^/^, but \i j3 it increases exponentially fast, being near unity when be^^/^^ ~ 1. 

In a general situation where the number of half-oscillations n is some function n{(f)Q), the change 
in the initial field necessary to change the number of half-oscillations by unity is 



dn((/)o) 



d(j)o 



(38) 



This needs to be compared with the width of the initial distribution cr ~ i^i in order to estimate the 
bias. If Hidn{(j)o) /dcpo ^ 1 the bias will be exponentially damped and the domain wall network will 
be stable. On the contrary, if Hidn{(f)Q) / dcpo <^ 1 the bias will be of order unity and the subsequent 



evolution will make the wall network collapse exponentially fast |16-18 



V. HISTORY OF THE FIELD ENERGY DENSITY AND WALL FORMATION 

The history of the energy density of the coherent field component for the Z2-symmetry breaking 
potential (|^) is approximately as follows. During inflation the energy density in the (j) held is of 
0{Hf), and the value of the field is 0o ~ 0.3A~^/^i?i in stochastic equilibrium. (We have assumed 
that the quartic term in the potential dominates so that Hi > h, where h is the height of the 
potential barrier; this is a necessary condition for inflation to produce domain walls, since otherwise 
the field will be settled in one of the two minima both during and after inflation.) If the symmetry 
is restored by thermal effects following reheating after inflation then domain walls will form again 
by the Kibble mechanism [Q] . However if the field is sufficiently weakly coupled this will not happen 
p8| and H\ < h will then be a sufficient condition for inflation to solve the wall problem]^ 

After inflation, the field perturbations do not evolve significantly until H becomes less than 
uj ~ ~ Q.'iX^/^Hi. Subsequently the field starts oscillating. At this point its energy density 

is still ~ QmHf while the energy density of the universe is ~ 3M| ~ O.SA^/^Ff Thus 
for the field (f) not to dominate the energy density we require A^/^Mp > 0.03/7f , i.e. 0o should 
not significantly exceed Mp, which seems natural. If A < \Q^^{H\/M-p)'^ the stationary stochastic 
distribution for the initial conditions will not apply. 

The energy density of the oscillating field redshifts as radiation hence its relative energy, com- 
pared with the total, either decreases or remains constant depending on whether the universe 
is matter-dominated or radiation-dominated. Therefore the energy density of the field can never 



^In principle symmetry restoration can also occur through non-thermal effects during "preheating" 
leading to the formation of topological defects [^] — however according to the results of Ref. |3^ 
this will not happen for the model under consideration here. 
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dominate since it did not do so initially when the field was began to oscillate at the end of inflation. 
However after the domain walls form, the relative energy density starts increasing again (at least 
until the wall network decays due to the bias). The field is released when t = ^ (A^/^ffi)~^ and 
subsequently its energy density decays as = {R^ / R)'^ . The walls will form when decreases 
below /i^ i.e. at a time determined by 

Using the results of the preceding sections we show in Fig. |^ the outcome for wall formation 
and survival in the parameter space of the scalar field model for two specific set of inflationary 
parameters — (a) Hi = 2 x 10^ GeV corresponding to an inflationary scale of ~ 10^^ GeV, 
Treh = 10^ GeV, and (b) Hi = 10^^ GeV corresponding to an inflationary scale of ~ 10^^ GeV, 
Treh = 10^ GeV. Reheating is assumed to occur while the inflaton oscillates in a quadratic potential, 
so the scale-factor evolves as for a matter-dominated universe and reheating ends at treh — -^p/^reh- 



VI. DISCUSSION 

We have developed the tools for computing the bias in domain wall formation in a specific 
Higgs-like model with a Z2-symmetry. As a general rule we find that the quantity which controls 
the bias is Hi {dn{4>o)/d(j)o), where the function n^cpo) gives the total number of half-oscillations 
performed by the field, starting with the value (po at the end of inflation, until it settles in one of 
its symmetry-breaking minima. This has to be evaluated taking into account both the form of the 
scalar potential and the expansion history of the universe. An exponentially small bias corresponds 
to Hi (dn(</>o)/d0o) ^ 1 while the bias increases approximately linearly for Hi (dn(</>o)/d(/>o) < 1. 
To obtain the probability distribution of the bias, one needs the initial probability distribution 



for (pQ. We have used the stochastic description |19| of the fleld fluctuations during inflation to 
compute this. 

The results of our detailed study for the Higgs-like potential are shown in Fig. |^ adopting both 
a high (10^^ GeV) and an intermediate (10^^ GeV) energy scale for inflation. We see that stable 
domain wall formation does not occur in most of the parameter space. For intermediate-scale 
inflation, the region where a problematic stable domain wall network forms is much smaller than 
for GUT-scale infiation and the domain wall problem is practically eliminated in this case. 

The results of this paper can be applied to other models as well. For a harmonic potential 
as in the case of an axion field we recover the results of Ref. 0]. For a periodic potential the 
distribution of initial values of d) is not relevant and we have a situation similar to the case of 



a Higgs-like potential in a radiation-dominated universe which was studied in Section IV. The 
parameter that controls the bias in this case is Hi/f where / is the period of the potential (i.e. the 
scale of Peccei-Quinn symmetry breaking for the axion field). 



In Ref. 1 12 1 the authors consider domain wall formation in a model with a Higgs-like potential 
for a dilatonic-type scalar field with a very small coupling constant, A ~ 10"*^^, and a vacuum 
expectation value v < 10^^ GeV. As the authors acknowledge, the issue of domain wall formation 
is somewhat subtle in this extremely weakly coupled theory. For the model considered in the 
present work, such couplings imply wall formation around the present epoch or later (see Fig.Q), 
so they would not in fact be astrophysically relevant. Furthermore the walls are formed in the high 
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bias region unless H[ ^ 10^^ GeV. However unlike the model considered in this work, the dilatonic 
field of Ref. |12| also couples universally with matter through the term 



>Cint = exp ( 0.^, (40) 

where 9^ is the trace of the energy-momentum tensor and M*^ ^ (10^ — 10^)Mp. During inflation 
the trace 9^ is non-zero, driving the field quickly to large negative values (p ~ —70AI* (where the 
exponential factor in Eq.(pOD makes the size of this interaction term comparable to the Higgs-like 
potential). Thus the effective mass is much smaller than Hi and the generation of fluctuations of 
size Hi is unavoidable. After inflation the term (|4^) decreases rapidly, so the field must relax in 
the quartic potential alone, starting at scales higher than the Planck mass. As we have previously 
discussed the absence of a force term strong enough to drive the field to the origin will make 
the potential energy of the field dominate the energy density of the universe. In Ref. the 
authors also introduce a term that couples the field coherently to the thermal bath during the 
radiation-dominated era, 

^therm = '^H^^l^^ ~ ^c/*' , (41) 

where k is a numerical constant. The intention in doing so is to drive the field to the origin, 
restoring the symmetry, so domain walls would be formed when H falls below the (vacuum) mass 
of the field. However, if cj) is initially greater than the Planck mass this term will exceed the energy 
density of the universe so the treatment is not consistent. This problem does not occur if k ^ 1 
but in that case the mass of the field is always much smaller than H so the thermal term cannot 
affect the field. 

Apart from this problem with the initial conditions, it is interesting to examine the effect of the 
coupling ( |4l| ) with the thermal bath on the damping of the oscillations. The equation of motion 
for the field in the radiation-dominated universe with such a potential term is exactly solvable, 
with the general solution 

For K ^ 1/4 the friction dominates and the field decreases very slowly as t"*^/^ (dominant term 
in Eq.(|4^)), while for k > 1/4 the behavior is oscillatory and the amplitude decreases as ~ 
R^^/'^ . Thus in either case the field amplitude decreases more slowly than for the quartic potential 
we have considered, where (j) ~ or even for the quadratic potential, where cj) ~ R^^^"^- This 

makes it even more unlikely that the walls can be formed before the present epoch. Thus this 
interesting attempt to do away with dark matter in galaxies by modifying the gravitational force 
law cannot work. 
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FIG. 4. The schematic situation concerning wall formation in the parameter space of the model 
dl]) discussed in this work. In region I the inflationary perturbations are too small to jump the 
potential barrier and the symmetry remains always broken. In regions II and III walls form with, 
respectively, high and low bias. In region IV the walls have not yet formed by the present time and 
the field is still oscillating. The walls form in the matter-dominated era above the upper dotted line, 
during the reheating epoch below the lower dotted line, and during the radiation-dominated era in 
the region in between. The left vertical dotted line separates regions where the oscillations begin 
at the reheating epoch (right side) or in the radiation-dominated era (left side). The right vertical 
dotted line separates regions where stochastic equilibrium applies (right side) or where different 
initial conditions apply (we have conservatively assumed (po ~ Mp for the initial condition in this 
region). The dashed lines indicate the domain wall tension a corresponding to wall domination at 
present {a^^^ ~ 100 MeV, lower line), and excessive CMB anisotropy {a^^^ ~ 1 MeV, upper line). 
Above these lines the tension is low enough to avoid conflict with observations. In the neighborhood 
of the line separating regions II and III, the bias takes values between zero and unity. In region 
III it approaches zero exponentially fast so the wall system survives. In region II the bias is large 
and the wall network collapses so there is no cosmological domain wall problem. The inflationary 
parameters are (a) Hi = 2 x 10^ GeV, T^eh = 10^ GeV, and (b) Hi = 10^^ GeV, Zeh = 10^ GeV. 
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